Casimir forces in a T operator approach 



00 

o 
o 

(N 
o 

Q 

(N 



> 
O 

o 
o 



X 



Oded Kenneth^ and Israel Klich^ 
(1) Department of Physies, Technion, Haifa 32000 Israel 
(2) Department of Physics, California Institute of Technology, MC 114-36 Pasadena, CA 91125 

We explore the scattering approach to Casimir forces. Our main tool is the description of Casimir 
energy in terms of transition operators. The approach is valid for the scalar fields as well as 
electromagnetic fields. We provide several equivalent derivations of the formula presented in Kenneth 
and Klich [Phys. Rev. Lett. 97, 160401(2006)]. We study the convergence properties of the formula 
and how to utilize it, together with scattering data to compute the force. Next, we discuss the 
form of the the formula in special cases such as the simplified form obtained when a single object 
is placed next to a mirror. We illustrate the approach by describing the force between scatterers in 
one dimension and three dimensions, where we obtain the interaction energy between two spherical 
bodies at all distances. We also consider the cases of scalar Casimir effect between spherical bodies 
with different radii as well as different dielectric functions. 
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I. INTRODUCTION 

The Casimir force [l[ is one of the fundamental pre- 
dictions of quantum physics. It explores the interplay 
between a quantum field and external " classical" like ob- 
jects such as boundary conditions, background dielectric 
bodies or space-time metric. While the classical objects 
modify the behavior of the field due to their presence, the 
field, in turn, acts on the objects, typically by exerting 
forces. Much work has been devoted to understanding 
the effect, as it appears in varied branches of physics: 
from condensed matter (interaction between surfaces in 
fluids) to gravitation and cosmology. 

The first precise measurement of the effect by Lam- 
oreaux , signaled a new age of Casimir force measure- 
ments, and led to a revived interest in the theory behind 
the effect. In recent years, the force between various ob- 
jects (such as two plates, plate and a sphere, corrugated 
plate and sphere, etc (0, 0, was measured. More- 
over the dependence on various properties of the mate- 
rials used , such as corrections due to finite conductiv- 
ity and temperature (@, @|)) as well as on geometry has 
been investigated. There is excellent agreement between 
the experiments and the theoretical predictions, which 
is being constantly improved. For an introduction to the 
subject as well as reviews of progress see e.g. 0, H B [l3 ■ 

Different variants (both material and geometric) of the 
force have been proposed, discussed and motivated by 
pure theoretical interest as well as by potential eventual 
application in nano-mechanical structures 
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The original method used by Casimir, that of mode 
summation, has led to a large body of work on the 
effect in simple geometries, where the modes may be 
exactly computed. For more general cases one has to 
use other available approaches such as the Green's func- 
tion approach or the path-integral approach. Signifi- 



cant progress in utilizin g th ese techniques numerically 



has been reported lately [13|, [14 1 



In the ID case, scattering approach to Casimir physics 
has proved very useful. Indeed, many of the calculations 
of Casimir interaction between bodies are based on scat- 
tering theory, as the photon spectrum in an open geome- 
try is continuous and it's description requires scattering. 

In this paper, we explore a scattering approach to 
Casimir effect in higher dimensions. The approach is 
based on analysis of a determinant formula for Casimir 
interactions obtained in Ref [lB|, and may be viewed as 
a generalization of previous formulas, especially related 
to scattering, such as the Lifshitz formula [lB|, and the 
results of Balian and Duplantier [l3|- Within this ap- 
proach, the Casimir energy is encoded in a determinant of 
the operator I — TaGqTbGo where Ta, Tb are Lippmann- 
Schwinger T operators associated with bodies A and B 
and Go is the photons Green's function; we shall there- 
fore refer to the formula as the TGTG formula. 

In [l^ it was shown how general results regarding the 
direction of the force between bodies related by reflection 
can be obtained from the TGTG formula. For example, 
the sign problem of interaction between two hemispheres 
was resolved. This result was subsequently extended to 
a large class of interactin g fie lds possessing the "reflec- 
tion positivity" property [l9| (See also [20|, where use 
is made of reflection positivity arguments to infer attrac- 
tion between vortices and anti- vortices in a frustrated XY 
model). In [21] an alternative derivation of the formula 
was presented. 

The paper is organized as follows. In section (II) we 
start with a derivation of the determinant formula, as 
well as supply alternative derivations in terms of Green's 
functions and the T operator of a pair of perturbations. 
Section (III) illustrates how one obtains the appropri- 
ate formula in the vector (electromagnetic (EM)) case. 
Sections (IV) and (V) cover simphfied cases: the special 
case of a body placed next to a perfect mirror, and the 
dilute limit, which deals with very weak dielectrics by 
expanding round e = 1. 
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We then proceed to show how the formula is to be 
applied in actual calculations. We explain how the for- 
mula is to be used together with partial wave expansions 
of the scattered states (sections (VI) and (VII)). In ID 
where only two modes (left and right movers) exist at 
each uj this leads to a known closed form formula for the 
Casimir energy in terms of reflection coefficients (see, e.g. 



2al23|) 



In Section (VIII) we use spherical waves to obtain an 
explicit expansion for the Casimir interaction between 
compact bodies. We demonstrate this by computing the 
force between two spheres at all distances, thereby gener- 
alizing the approach of 2^ to spheres beyond Dirichlet 
boundary conditions, and going beyond the proximity 
force approximation. We also consider cases of spheres 
with unequal radii, as well as as spheres with arbitrary di- 
electric function. In section (IX) results are extended to 
the electromagnetic case. While this manuscript was fi- 
nalized we have learned that related results were reported 
in [2l| for the Casimir interaction between spheres with 
equal radii, as well as an alternative derivation of the 
determinant formula. Our results perfectly agree with 



those of [2J] and [21|. 

A number of appendices describe some technical details 
of the calculations. Most notably, in appendix (B) we 
give additional details about the mathematical validity 
of the formula, which were not included in [l^. These 
details, help establish for the first time rigorously the 
validity of the present approach to calculations of Casimir 
forces. In particular we show that the formula is given 
in terms of logdet(l -I- A), where Tr|^| < oo, and so 
mathematically well defined. (This appendix is written 
in a "mathematical physics" style, and may be skipped 
by readers not interested in these issues) 



II. THE TGTG FORMULA: CASIMIR 
INTERACTION AS A REGULAR 
DETERMINANT. 

In this section, we explain how the part of the free en- 
ergy of a Gaussian theory that depends on distance be- 
tween bodies, and as such is responsible for the Casimir 
force, may be expressed in terms of a regular determi- 
nant, and discuss some of its properties. Some of the 
material covered here appeared in the literature, how- 
ever, as far as we know, the final formula was never writ- 
ten in this general form; furthermore, it's mathematical 
properties where not rigorously addressed previously. We 
note, however that an elaborate and rigorous analysis of a 
related problem involving impenetrable discs was carried 
out in [25|. 

We start by presenting the derivation of the de- 
terminant formula in the path integral approach 
2^ , 26l . 27 1 . We first treat the case of a scalar field and 
explain later how the result is extended to the EM field. 



Alternative derivations of Eq ((TT|) are elaborated in the 
following subsections. 

The action of a real massless scalar field in the presence 
of dielectrics can be written as 

S[^] -Ud'r f ^0:(V2 +c.2,(x,c.))0. (1) 



where 6*. = 



and e{uj, x) = 1 + x(x, lu) is the dielec- 



tric function (we use units h = c = 1). This action is 
the simplest action which yields the scalar analog of the 
Maxwell equation 



V X V X A -e(uj,x)A 







(2) 



for the vector potential in the radiation gauge. Alterna- 
tively this action can be derived by coupling a scalar field 
to an auxiliary field living on the regions of space where 
e ^ 1, and then integrating out these fields, as done, e.g. 
in [27]. 

Formally, the free energy of the system is obtained from 
the partition function Z given by: 



V 



js\< 



(3) 



Performing the Gaussian integration, one finds that the 
change in energy due to introduction of x in the system 
is 

Ec = E^- E^=o = (4) 
-^ ^ log detA (1 + c.2x(x, Lu) (V^ + ,.2 ^0)-^) . 

At this point, we encounter one of the main features of 
Casimir physics - the need to properly isolate the phys- 
ically relevant part of the energy out of a formally ill 
defined expression. A determinant (such as in Eq (jl])) 
is mathematically well defined only if it has the form 
det(l + A), where A is a "trace class" operator, i.e. 
^j|Ai| < CO with Xi eigenvalues of A (For properties 
see appendix [Cl and ISj, |28|). If A is not a trace class op- 
erator, one may obtain different or infinite results for the 
determinant, depending on the order in which the eigen- 
values of 1 + A are multiplied. The expression above is 
not of the required form. To sec this note that 'A' in our 
case is given by: 



w\(x,w)(V2 + +iO)" 



(5) 



This is an operator of the form g{x)f{V). If such an 
operator is "trace class" then it's trace is known to be 
given by the Birman-Solomyak result [28l|: 



Tr(g(a:)/(*V)) = / d'xgix) / d'^kfik) 



(6) 



in our case we have J d^xxix) < oo, however J d^k{—k^+ 
-|- iO)~^ diverges, and is indicating that the operator 
involved doesn't have a well defined trace 1361. 



3 



A 





FIG. 1: Bodies A and B 

As such, the expression (j4]) only has meaning when 
specifying physical cutoffs. Removing physical cutoffs 
will leave us with an ill defined expression and so we 
keep in mind cutoffs at high momenta in the notation 
detA. 

At high frequencies, x(ct;,x) — > provides a physical 
frequency cutoff. For Re w, Im w > both x(^) ^-nd 



(V2 



iO) ^ are analytic, justifying Wick-rotation of 



the integration to the imaginary axis iuj ending up with: 



c 



AljJ 

2^ 



logdetA(l + w'x(x,iw)Go(x,x')) (7) 



Where G'o(x, x') = (x| |x'). Restricting the op- 

erator (1 + Lo'^xGo) to the support of x (more precisely 
to L'^{Supp{x))) clearly does not affect its determinant. 
Note that Eq. ^ is still ill defined if one removes the 
cutoff, as can be immediately seen from the argument 
based on Eq. ©. 

We now consider the case depicted in Fig. [1] of two 
bodies A, B immersed in vacuum, x is assumed nonzero 
only inside the volumes of the two dielectrics A, B and 
we therefore consider in the following (I + w^xCo) as an 
operator on Ha ® Hb Ha ® Hb where Ha = L"^ {A) 
and Hb = L'^{B) js^]. It is then convenient to write 



(I+w^xGo) 

1a + u?xaGqaa w^xaGoab 
^^XbGoba ^b+(^^XbGobb 



(8) 



where Gqq^ is Go restricted to Ha Hp (equivalently, 
Goa/3 = PaGoPp, where = I e and Pg = ® I, are 
projections on Ha,Hb respectively). It turns out that 
the part of the energy that depends on mutual position 
of the bodies, and as such is responsible for the force, 
is a well defined quantity, which is independent of the 
cutoffs. To see this, we subtract contributions which do 
not depend on relative positions of the bodies A, B: 



Ec = Ec{A\Jb) - Ec{A) - Ec{B) 



(9) 



As in Ref. 27], this amounts to subtracting the diagonal 
contributions to the determinant, which are not sensitive 
to the distance between the bodies, (i.e. only contributes 



to their self energies). This yields 



Ec — 



(10) 



/■oo dcj / i„„ I^+^'^XaGoaa ^'^XaGoab 

Jo 27liogdetA(^ u;'xbG,ba 1+u;'xbGobb, 

, , ^ / 1 + uj^xaGoaa \ \ 

log detA ( , - ■^--^Go^JI 



1 + lu^xbG, 

/ri^(iogdetAf J, ^; 



1 TaGoab 
b^qba 1 



where = 



I X 



det 



Y \ 



Xa- Finally, using the relation 
= det(l - FAT), 



which holds for block matrices we have: 



Ec{a) 



rf^logdet(I 



TaGoabTbGoba)- (11) 



Henceforth, we refer to pT|) as the TGTG formula 
throughout the paper. Up to Wick rotation, the oper- 
ators Ta are exactly the T operators appearing in the 
Lippmann-Schwinger equation, as will be discussed in the 
next subsections. The Wick rotation T{llj) T{iuj) has 
the effect of turning Tq, into hermitian operators as well 
as of avoiding potential singularities (which may occur 
at real frequencies). 

In Eq. ([TT|) . we disposed of the cutoff A as the ex- 
pression is well defined in the continuum limit. In prac- 
tical terms this means that replacing the infinite dimen- 
sional matrix of I — TaGqabTbGqba by its its upper- 
left n X n block with n large enough and calculating 
the resulting ordinary determinant, gives an arbitrarily 
good approximation to a (finite) quantity which we call 
det(l — TaGoabTbGqba)- This point is discussed in de- 
tail in appendix [Bl where we prove some mathematical 
properties of the operators involved. The details are not 
essential for understanding the applications of the for- 
mula, so a reader not interested in mathematical rigor 
may skip them. 



Dirichlet and Neumann boundary conditions 

In many cases, and indeed in the original presentation 
by Casimir, one is interested in sharp boundary condi- 
tions, such as Dirichlet or Neumann. Sharp boundary 
conditions result in singular energy density at the surface, 
as field modes are required to vanish for all momentum 
scales. Typically, the local energy density diverges as the 
inverse fourth power of the distance from the boundary 

It is important to point out that the above consid- 
erations also describe the conducting case with minor 
changes. Following [2^, assume conducting boundary 
conditions are given over a surface S, parameterized by 
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internal coordinate u and by the embedding in R'^ given 
by x(u). One may describe a simple metal by taking 

x(ia;) = ^^^2 on S, and letting E have a thickness of a 
few skin depths l/ujp, I ^ 0(1), here LOp is the plasma fre- 
quency (proportional to the effective electron density in 
the metal). In the limit of large LOp one retains the same 
expression as ()lip . with the following substitutions: 

Ec{a) = (12) 
2^ /o°° dt^(logdet(l - MsAT+^A^AST+^) 
where in the Dirichlet case M. is given by: 

M(^)(u,?/';t^) = Hy5RGo(x(7/),x(u'))v/5R) (13) 

and acting on the surfaces S. Similarly Neumann bound- 
ary conditions may be treated in the path integral 
method by taking |30 |: 

mW(m,u';w)= (14) 
\J 9{u)^ g{u')dn(^u)dn{u')GQ{y^(vL),yi{u')) 

in We remark, that rigorous discussion of the for- 

mula in the Neumann case requires further analysis which 
we did not pursue in this paper (see remarks after Eq. 

Derivation using Green's functions and T operators 

To make contact with Green's function approach we 
supply in this section an alternative derivation of the 
TGTG formula. Most of the derivation is standard and 
may be skipped by readers interested only in new results. 
However, we point out that our approach where the T 
operator of combined scatterers are utilized seems new. 
Here, we use the Green's function in order to express the 
density of states (DOS) of a differential operator with 
background, and then perform the mode summation by 
integrating over energies. 

We briefly remind the reader some of the required ma- 
terial. The standard discussion of this is usually done in 
the context of non-relativistic quantum mechanics. The 
retarded/advanced are then defined by: 

(E±is-n)G^{E)^I (15) 

This equation should be understood as an operator iden- 
tity. If 7i is a differential operator, for example 7i = — A 
then it is the operator form of the differential equation: 

(E±is + /\)G{x,x')^8{x-x') (16) 

Using the representation (n|G*(£')|n') — 
lim„^n 1 one then finds that the DOS is 

^ ^ tj±%S— tjn 

given by: 

ilmTrG±(i?) = ^ V <5(i? - i?„) - TpiE). (17) 

TT 

n 



Noting the identity 

dE log[i? ± - = (18) 
E ±is — En 

one can rewrite this as 

p{E) = ±-lmdETr\ogG'^{E). (19) 

TT 

We are more interested in the relativistic version of 
this. (Indeed the Casimir force vanishes in the non- 
relativistic limit, as the exchange of very massive virtual 
particles is suppressed.) In the relativistic context the 
Feynman propagator G is defined by a similar formula to 
that of G+: 

H{lj^ +is)G ^ I (20) 

For example, the action ([T]) corresponds to H — — A — 
Lu'^e. In free space e = 1 we obtain the same equation as 
Eq. (|16p apart from the substitution E ^ uf^ . (There 
is also a not very interesting conventional overall minus 
sign, which is the reason some signs in the following equa- 
tions are different from what the reader may remember.) 
In the presence of nontrivial background (e.g. dielectric) 
the uj dependence of H can take quite an arbitrary form, 
a fact that slightly complicates the derivation of the DOS. 
One may take advantage of the relation 

where F{x) is any real function having simple zeroes 
at the points {xn}. Indeed, away from the zeroes 
the fact that F is real guarantees vanishing of 

the l.h.s while near the zero Xn we have Im-^r^^^^ = 

I"i(J±lSffc) = ^ T^'5(x - x„). Gen- 

eralizing the relation from real functions F{x) to hermi- 
tian operators H{uP') [s^ allows writing 

Ima^TrlogG(w) = -ImTriJ'(w)G(w) = Tip{uj) (21) 

which is the obvious analog of Eq. (fT9| . (Note however 
that similar generalization of Eq. (|17p would usually be 
false.) In ((2T|) we implicitly assumed a; > to avoid an 
extra sign{uj) factor. 

Now, assume that Go is known for Ho and we add a 
perturbation V , i.e. 

{Ho{lj^ + i3) + V{uj'^ +is))G = L (22) 

The change in DOS due to introduction of the potential 

V is formally: 

Ap = -Im 9^ Tr log GGo^ (23) 

TT 

We will be interested in the change in energy due to 
change in the distance a between two separated po- 
tentials Va and Vb, which make up V. So we take 

V = Va + Vb. 
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Thus, 

da^p{Lo) = ilm5^a„Trlog(GGo-^) = (24) 
ilnia^9,Trlog(/ + Goy)-i. 

Defining the T matrix by 

T =V{I + GoV)-\ (25) 

we may also write 

da^p{Lo) = ilm9^9,Trlog(/ - G^T). (26) 

TT 

Alternatively, formally writing "9adet(Vyi + Vb) = 0" 
since Va , Vb act in different subspaces one can write that 

daApiij) ^-Imd^daTilogT. (27) 

TT 

this last formal expression, however, should be handled 
with care, and so we avoid using it. 
The T matrix satisfies: 

G{uj) = Go(c^) - Go{uj)T{uj)Go{u;) (28) 

(here cj > is actually lo + is), and frequently appears 
in scattering theory. Also note T = V — VGV. 

The operator T appears in the Lippmann-Schwinger 
equation as follows. Given a solution (j) of the free equa- 
tion, without a potential Ho{uj)(p = 0, one constructs a 
solution ip of the eigenvalue equation {Ho + V)ip = hav- 
ing the same incoming part V'm = (f'in- Formally, this is 
done by looking for a solution of: 

^ = GqV^P, 

which is the Lippmann-Schwinger equation. It follows 
that ip = {I + GoV)-^(j) = (I-GoT)^, thus we may build 
a new solution from a solution </> of the free equation. 
For example, choosing (/) to be a plane wave solution, one 
obtains 

V'fe = 6*'==" - J dk'Goik'){k'\T\k)e''''''. (29) 

Note that our relativistic normalization convention 
implies that T is related to the scattering matrix via 
S =1- 2TTiS{uj^ - Ho)T. 

We now address the case of two potentials Va,Vb- We 
assume for simplicity that cutoffs are in place, and so 
work with the T operators as matrices. We compute 
the joint transition matrix for both perturbations Ta\jb, 
and show that the part independent on "self energy" is 
exactly Eq. (fTTj) . 

By using the formula (|28l) as Gi = Gq — GoTiGo (with 
i = A,B), together with the definition of T (|25|) . and 
straightforward algebraic manipulations we obtain: 

l + Go{VA+Vn) = (1 ^ ^oTa) i_^GoTbG„Ta ~ GoTbGSO) 



and so the joint T operator of a pair of perturbations 
may be factored as 

Taub ^ (Va + V^i3) i+G„(y^+y^) = (31) 

{Va + Vb){1 - GoTA) ,_aJ,GoTA ^ ~ ^oTs). 

The important feature of this expression is the observa- 
tion that the only part of the expression which directly 
mixes between the A and B is the factor 1 — GqTbGqTa- 
Indeed, plugging Eq. dSO]) in Eq. (l24|) we see that the 
contribution of frequency lu to the force is now given by: 

daAp{Lo) = ilm9^9,Trlog(/ + GoiVA + Vb))-' - 
iW^aa[logdet(l - GoTa) + logdet(l - GoTb) - 
\ogdet{l~GoTAGoTB)] = 
-^Imd^da logdet(l - GoTaGoTb) (32) 

leading again to our expression for the energy Eq. (jlip . 

Alternatively one may simply verify the correctness of 
Eq. PT|) by noting that: 

1 - GqTaGoTb = 1 - GoVa 1+gVa ^"^B 1+GoVb = 

{(1 + GoVa)(1 + GoVb) - GoVaGoVb} 
^ i+GoVs = i+g„Va [1 + '^oIVa + Vb)] i+g\vb 
and using Eq. 

III. THE ELECTROMAGNETIC CASE 

Here, we follow the approach of [sij. The statistical 
properties of the electromagnetic field in a medium are 
described by the appropriate photonic Green's function. 
The electromagnetic fields are derived from the electro- 
magnetic potentials A", a = 0, ..,3. (It is convenient 
to work in the gauge A'^ — 0.) The retarded Green's 
function Vik is defined by: 

I?,fc(Ai,A2)= (33) 

r {A,{X,)Ak{X2) - Ak{X2)A{Xi)) h < h 
I otherwise 

where Xi,X2 are 4-vectors X = {X^,..X^) and k,i — 
1,..3. The angular brackets denote averaging with re- 
spect to the Gibbs distribution. 

The interaction of the electromagnetic field with a clas- 
sical current J put in the medium is given by 

Kubo's formula allows us to treat this interaction 
within linear response. By Kubo's formula the mean 
value Ai in presence of a current J satisfies: 

A:(r)^ ^-^cj 2?^'fe('^;r,r')Jfc(r').d3r', (34) 
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where 



Vg{t;r,r')dt 



(35) 



The function V is sometimes referred to as the general- 
ized susceptibility of the system [sij . 

From Maxwell's equations it follows that in a medium 
with a given permittivity tensor 6^ , permeability tensor 
fj,ij, and current J, the vector potential Ai satisfies: 



(V X (m^Vx) - ^e)A = —J 



(36) 



Substituting Eq. (134)) in Eq. ([36)) . we see that I? is a 
Green's function for the equation: 



V X /i-iy X P - —eV = -47rM5(r - r') 



(37) 



where I is the three dimensional unit matrix. In the 
following we shall work in units where c = h = 1. 

The Green s function V is then used to obtain the well 
known expression Eq. (80.8) in Lifshitz and Pitaevskii 
3l| for the change in free energy due to variation of the 
dielectric function e at a temperature T: 



SF 



1 



OO 



SFo + -T ^ 



(38) 



n— — OO 



Here is the free energy due to material properties 
not related to long wavelength photon field, and a;„ — 
2TrnT are Matsubara frequencies. V is the temperature 
Green's function of the long wave photon field given by 

Eg ((38)) may be written as SF = SFa + SFc where 
Fc = 5T^Er=-oo[logdetA(V X V X +c^2e(a;, zc^„)) 



logdetA(V X V 



;)] 



|T Er=-oo logdetA(l + ujIx{x, iuJn)T>o{iuJn))- 



Here 2?o(2^, x\ iLUn)i 



Note that Fr 



is exactly the same as ([7]), with the scalar propagator Go 
replaced by the vector propagator Vq. For later reference 
we write here the explicit expression for Vq: 



T)o,j (fc, iu) 



An 



fc2 



r(4 



(39) 



Thus, starting with this expression, one repeats Eq. ^ 
and Eq. dTU]) to get Eq. replacing Go by Vq ev- 

erywhere (including in the definition of the T operators). 
The analysis of the determinant now proceeds exactly as 
in the scalar case. 

Alternatively, the EM case may similarly be derived 
starting from the functional determinant corresponding 
to an EM action analogous to Eq. ([T]) . In the axial gauge 
Aq = this action takes the form: 

S^\U^v I X V X We(x,u;))X. (40) 



A permeable body may similarly be described within 
our approach by replacing the dielectric interaction term 
Aio'^xix)^ in the Lagrangian by a magnetic term: x 
(1 — i)V X One may then go on through our deriva- 
tion using the (differential) operator V x (1 — ^j^)Vx 
instead of —Lo^xi.^) everywhere. One major difference 
between the two cases is worth noting: whereas the di- 
electric term is always described (after Wick rotation) by 
a positive operator, the operator in the magnetic term 
turns out to be negative (for /i > 1). This fact can be re- 
lated to the known Casimir electric-magnetic repulsion. 
Moreover, the ideal n oo limit is seen to correspond to 
a Lagrangian in which the term x A)'^ is missing (in- 
side the body) which makes a highly irregular lagrangian. 
Analogy with a scalar field satisfying Neumann bound- 
ary conditions suggests that this situation may be de- 
scribed by dropping the (Vi/))^ term inside the Neumann 
body. There are also other arguments in favor of that 
approach [s^l, however we did not bring these argument 
to a completely rigorous form. 

IV. DIELECTRIC IN FRONT OF A MIRROR 

A somewhat simplified, but useful in practice, version 
of our formula is obtained in the case of a body placed 
close to a mirror. Consider the body A to the left of a 
Dirichlct mirror B located ed, Xn ^ cl/2. It is well known 
(using the image method) that the effect of the Dirichlet 
mirror is to replace the free propagator Gq by 



Gb(x, x') = Go(x, x') - Go(x, J(x')) 



(41) 



where J(xj| ,x±) = (x|| , a—x±) denotes reflection through 
the mirror plane. This may be written as Gb — Go = 
—GqJ where J is the operator defined by ^^^{x) — 
ip{J{x.)). Noting the standard relation (1^5)) Gb — 
Gq — GoTbGo between the Green function in the presence 
of scatterer B to its T matrix one concludes GqTbGq — 
GqJ which when substituted in pip gives 

/>00 1 

Ec{a)^ -^^logAetil ^GoJTa). (42) 

An alternative (though closely related) approach is to 
note that by complete analogy to Eq. ([7]) the energy it 
costs to place a body A near a mirror B is 

Ec ^ I -^logdetA(l +w2xA(x,icj)Gs(x,x')). 
Jo 

Subtracting the energy Ec ~ ^ logdetA(l + 
Ci;2x^(x, iCi;)Go(x, x')) it cost to put A in vacuum then 
gives the Casimir interaction energy. Using the relation 



{1 + GbVa)/{1 
1 + (Gb - Go)T^ 



G0K4) = 

: 1 - G(,JTa 



(43) 
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leads again to ([42]) . 

Yet another way of obtaining the same result is by 
substituting xb — XS{xn — a/2) in the definition of Tb 
and doing the algebra. One then finds 



{x~x')j 



GqTbGo — 

^-q\x„\~-q\x'J 



(44) 



2q(Xu'-^+2q) 



which in the limit A — *■ c» reduces, as expected, to the 
expression GqJ obtained through the image method. 

We now address the case of a Neumann mirror. Note, 
that the Green function in the presence of a Neumann 
mirror is G = Go + Gt^J . By repeating the arguments 
above we find that the Casimir interaction between an 
object A and a Neumann mirror is given by a similar for- 
mula to involving the determinant det(l + GoJ'Ta). 
We remark, that while the Dirichlet mirror may be 
considered as the limit A — > cxd of a dielectric having 
e.g. Xb — A^(a;„ — a/2) (or in more realistic model 
Xb = ^0{xn — a/2)) it is hard to find a simple analog 
XBix) that would lead in a similar limit to a Neumann 
mirror. (See however remark at the end of the previous 
section). 

A similar treatment is applicable in the more physi- 
cally relevant EM case. The boundary conditions = 
may be enforced by requiring the vector potential to sat- 
isfy J^A = —A where is defined to act on vectors as 
JA{x) = {A\\{J{x)),-Aj_{J{x))) (Here A\i,A_l denote 
the components of A parallel and normal to the mirror 
surface. The temporal component is considered as a par- 
allel component though in practice we usually choose a 
gauge where it vanishes.) 

The EM Casimir interaction between a dielectric and a 
mirror is then given by a formula similar to Eq. (I42p 
with Go , replaced by the EM propagator Vq and the 
vectorial J defined above. 

It is interesting to also consider an ideal permeable 
mirror (having /i — s- oo,e = 1). This corresponds to the 
boundary condition = which may be enforced by 
requiring the vector potential to satisfy J^A = +A. Thus, 
the Casimir interaction of body A with such a mirror 
will be given by an expression involving the determinant 

det(i + rAi?o:^). 



V. DILUTE LIMIT 

In the following sections we consider strategies of using 
the TGTG formula in actual calculations. A particularly 
simple case is when x is small, which is commonly re- 
ferred to as the "dilute" case (and sometimes as "low 
contrast"). Here we briefly sketch how to best use the 
formula in this limit. As shown in the appendix (Theo- 
rem |BT61), one always have \ \^/GTGT^/G\ \ < 1, therefore 
we may expand the logdet(l — ..) expression pip in pow- 



J 2n ^-^ m 



(45) 



In the dilute limit Xa < 1, so one may also substitute the 
expansion 



T^ = -Y,i-^^XcGoru;\a 



(46) 



n=0 



in Eq. (j45p and compute the involved integrals to desired 
order. This expansion is the continuous equivalent to 
summation of two body forces, and is equal to the Born 
series appearing, for example, in Ref. 3l( . 



VI. SCATTERING APPROACH 

As remarked above, the operator TaGqTbGq appear- 
ing in our formula is closely related to scattering data. 
The purpose of this section is to clarify this relation and 
make it more explicit. In order to keep better touch with 
conventions used in scattering theory, we usually avoid 
in the following sections using Wick rotation and thus we 
work in Lorentzian rather then Euclidean space with real 
rather than imaginary frequency and with the Feynman 
rather than the Euclidean propagator. 

As mentioned above, the arguments of Go in Eq. (fTT|) 
never coincide, implying that when GQ{xa,Xb) is con- 
sidered as a function of Xb alone it is a solution of the 
(homogeneous) free wave equation. Thus one may ex- 
pand Goixa,Xh) in the form J2Caf34>ai^a)4>i3ixb) where 
{4'a{xa)}, {(t>f3{xb)} are some sets of free wave solutions 
of energy uj. There is, of course, great freedom in choos- 
ing the sets {(j)aixa)}i {(f'pixb)}- In practice one would 
choose these in a way that makes subsequent calculations 
easier. Since we consider TaGqTbGq as acting only on 
the volume of object A, these considerations also apply 
to the propagator on the right of this expression. 

The Lippmann-Schwinger operator T[uj) is related to 
the S-matrix by 



1 - 2Tri5{uj'^ - uj'^)T^ 



(47) 



Therefore, T(uj) has the property that its matrix element 
{a\T\[3) between a pair of free states a^(3 having energy 
Lo is equal to the corresponding matrix element of the 
transition matrix. Since the operator Tb in TaGqTbGq is 
sandwiched between a pair of free Feynman propagators 
corresponding to energy uj, we may identify it with the 
corresponding transition matrix. Due to the cyclicity of 
the determinant det(l — TaGqTbGq) the same is true of 
Ta. 

Substituting the expansion GQ{xa,Xb) — 
Y^Ca(3(t)*a{xa)(f>f3ixb) wc arrive at 



TaGqTbGq — 



E 



TA\a)C^p{(3\TB\P')C^'^'{a' 
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The Casimir interaction will then be given explicitly by 

duj 



E : 



27r 



logdet(l - K{iuj)). 



(48) 



Here Ka"a' = {TA)a"aCap{TB)p[3'Ca'P'- 



VII. PARTIAL WAVES EXPANSION 

In the following sections, we consider strategies of us- 
ing the representation (j48|) by restricting the K matrix 
to a finite subspace which gives the dominant contribu- 
tion to the force. Indeed, in many cases of interest only 
a few partial waves are significantly scattered; the best 
example for this is when objects are far apart, and from 
a large distance look point like. At this limit one expects 
significant contribution only from s-wave scattering. In 
the more general case, K may be approximated by a fi- 
nite dimensional matrix corresponding to several partial 
waves. In order to see how this works in practice we 
consider below a few simple cases. 



One dimensional systems 

A particularly simple case occurs when the system is 
one-dimensional. Consider, e.g., a scalar field in ID. All 
states of energy lu are then spanned by two modes: left 
and right movers — -^e**"^. Hence, in this 

case, the determinant Eq (fTTj) can easily be evaluated. To 
see how this is done, we write the Feynman propagator 
explicitly as 



Go 



dk 



27r Cj2 - fc2 + iQ 



2uj 



(49) 



We consider a pair of scatterers A, B such that A is on 
the left of B. This immediately implies that we have 
Xa < Xb and therefore 



GoBAixb, Xa) 



Juj{x^-Xa) 



2UJ 



-2iTT 

2lo 



\R){R\ (50) 



Similarly, we also have Gqab — Using this 

we see that the operator K in Eq. turns into the 

c-number 

K - {^)'{R\Ta\L){L\Ts\R) - fA{i^)rB{io). (51) 

Here tb (m) is the reflection coefficient for a wave hitting 
scatterer B from the left {A from the right) to be reflected 
back. Note that the normalization of T implied by Eq (l47l) 
is responsible to the cancelling of the factor ^^j^ • (Had 
we used relativistic normalization for |L, R) the factor 2lu 
would not have appeared.) We thus conclude 

det(l - TaGoTbGo) = 1 - rAiLo)rBiiu). 



The tilde on va serves to remind us that it is the reflection 
coefficient from the right side of A. 

We remark that f^(o;)rB(w) implicitly depends on the 
distance between A, B through the (phase) dependence 
of rA,rB on the scatterers locations. To make this ex- 
plicit, note that moving a scatterer a distance a affects 
the reflection coefficients as r — e~2""r, r — > e^^^f . 

Moving the scatterers a distance a apart therefore re- 
sults in 

det(l - TaGoTbGo) ^ (1 - e2--f^(c^)rs(c^)). 



Substituting in (I48|) we obtain the familiar formula for 
Id Casimir interaction between scatterers 22. 23l 33j. 



Multi-component field in Id 

The considerations used above for a single scalar fleld 
in one dimension extend to a situation where (j) — 
(01,02, ■■■4>n) is an n component fleld. In this case, the 
reflection coefficients ta^b turn into n x n matrices and 
one finds det(l - TaGqTbGq) = det(l - fA(w)rij(w)) 
where the determinant on the right is of a usual n x n 
matrix. 



Plane wave expansion. 

In physical three dimensional space, there are many 
different possible ways to expand the propagator 
Go{xa,Xb) = J2^at3<t>a{xa)4>/3{xb) in tcrms of free wave 
solutions {(j)a{xa)}, {4>/3{xb)}. In the next section we de- 
scribe the expansion in spherical waves (which is proba- 
bly the most useful expansion), and we demonstrate its 
use to calculating the Casimir force between compact ob- 
ject. However, for the sake of simplicity we first describe 
here a plane wave expansion which is the immediate gen- 
eralization of Eq. (j50p . A simple heuristic way to arrive 
at this generalization is to formally think of the field 4> 
in three dimensions as one dimensional field having in- 
finitely many components labelled by its transverse mo- 
menta. Indeed, such point of view has been successfully 
used in describing transport in quasi ID conductors in 
mesoscopic physics, whereby each transverse component 
corresponds to a scattering channel (see for example Ref. 
[3^). This suggests splitting k into its z-component k^ 
and its transverse components fcy — [kx,ky). The 3d 
propagator may then be written as: 



Go 



(2n) 



2ky^ 



Here 



' — fc| + zO may be either real and positive (for 

Lo^ > fcy ) or pure imaginary (for < fc^) in which case 
the iO prescription implies that it must be chosen on the 
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positive imaginary axis. Assuming that A is located to 
the left of B along the z-axis it follows that 

TaGoTbGo = (52) 
^{{qx,qy,-qz)\TB\{kx,ky,k^)):^{{k.^,ky,k^)\, 



where 



qz = 



qx — qy + *0 and k^ 



■-kl~k^ + lO. 

When considering only the terms satisfying uj'^ > 
ql + qy,kl + ky, Eq. ([5^ indeed looks hke a straight- 
forward generalization of the Id result. However, as this 
expression shows, to get the correct result one must also 
include the contribution of evanescent waves (qjj > uj'^). 
Upon Wick rotation, however, the distinction between 
ordinary and evanescent waves disappears. It may also 
be noted that (since in general qz ^ kz) the variation of 
the {{qx, qy, -qz)\TB\{k3:, ky, kz)) matrix elements upon 
moving B along the z-axis is considerably more compli- 
cated then in the Id case. 

The above representation may be helpful in problems 
where the scatterers A, B have exact or approximate pla- 
nar geometry (e.g. corrugated plates). Though the the- 
orem guaranteeing finite trace does not apply for infinite 
plates, one may show that dividing by the plate area 
leads to finite result. We remark that actual calculation 
of the determinant requires discretizing fc|[ which corre- 
sponds to assuming large but finite plates. Alternatively 
one may use Eq. (j45p with continuous fcii . 



VIII. SPHERICAL WAVES EXPANSION 

When describing interaction between two compact 
bodies, often it is convenient to represent the transition 
matrices T in a spherical wave basis. To do so, we choose 
two points Pa,Pb inside bodies A,B respectively. We 
parameterize the points of body A by the radius vector 
f — fX measured from the point Pa and the points of 
B by the radius vector r' = fs measured from the point 
Pb ■ The vector connecting Pa and Pb will be denoted 
by a (FigUl). In the scalar case, the free spherical waves 
centered at Pa , Pb are given by 



\{lm)A,B) = 



■jl{t^rA^B)Yi„i{fA,B), 



(53) 



with the normalization {u'l'm'\ulm) = 5mm' 5{lo — lo') ■ 
To use Eq. (^5)) . the scalar 3d Green function Gq = 
— is expanded in terms of the spherical harmonic 
functions centered around Pa and those centered around 
Pb- 



E 

Imd'm' 



\{lm)B)Ci {{l'm')A\ 



(54) 




FIG. 2: Coordinate system used for the partial wave approach 

where (See appendix El for a proof of the following equa- 
tions.) 



2ui 2-1 1" .m" ^ 



I I' I" 



m m m 



{co) 



(55) 



z'"+''-'4V(c.a)yr'™"(a), 



Yim are spherical harmonics, ji,hi are spherical Bessel 
and Hankel functions, and the coefficients 



C 



I I' I" 
m m' m" 



have known expressions in terms of the 3j symbol or as 
an integral of spherical functions: 



C 



I V I" 



m m m 



= (-1)™ ^47r(2/ + l)(2;' + l)(2/" + l) X 

I I' r\ { I v I" 







(56) 



C;m;/'m'(«w) = — gimd 



In actual computations, it is often more convenient to use 
the Wick-rotated expression. This may be expressed as 

Where the coefficients 

(57) 

'-;^Ki„+l{uja)Yi„ra"{a), 



are real. Here Kl,,_^_l are modified Bessel functions of the 
second kind. Equations (|55l57p may be somewhat simpli- 
fied by choosing the z-axis along a. 

The above expansion of G^^ allows expressing 
TaGqTbGo in terms of matrix elements {l'm'\T\lm) of 
the transition matrices of the two scatterers. The Casimir 
interaction may then be written as in Eq. (j48p where 

Klrn-l'm' = (58) 
( 1) ^'^ " (^A)/m;/i mi^/i mi ;/2*Ti2 (-^-B )/2^^2 ;^3^3^^3^^3 ■ 



Here Cim-i'm' is given by Eq. fMl) or Eq. ([57]), 
summation over /i, toi, ^2, ^3, ^3 is implied and we 
note that the extra sign resulted from C/m;;'m'(— a) = 

(—1)'+' Clmd'm'ia) = Cl'm';lm{o-)- 

If we assume that only waves having I < Iq are 
significantly scattered then K will turn into a finite 
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{Iq + 1)^ X {Iq + 1)^ matrix (since the dimension of the where we used the notation: 



subspace I < Iq is J2iLoi'^^ + 1) = ('o + 1)^)- We stress 
that this argument does not require us to assume spher- 
ical symmetry of the scatterers. 

When A, B are very far apart, the interaction between 
them is governed by waves of very low frequency and 
therefore also low I. At this limit the leading contribu- 
tion comes from the s-wave scattering transition matrix 
element {I — 0\Ta,b\1 = 0) — '1ijP'\a,bI'^ ^ where A is the 
scattering length. 

The matrix K then reduces to the scalar K = 

-w^a^Ab (/io^^(^«)) = 47r^^^e2*'^'^. Doing the in- 
tegral ([i5|) one arrives at 



Ec = 



This limit corresponds to the scalar version of the well 
known Casimir-Poldcr interaction. Our formalism, how- 
ever, allows calculating corrections to it up to any de- 
sirable finite order in i. For example for two Dirichlet 
spheres of radii i?i , i?2 at distance a between their centers 
the expansion gives: 



E = 



R1R2 _ RiR.2{B.i+R2) 
-Ri R2 (34_R^ +9fli -R2 -l-34_R2 ) 
-Rl-R2(-Rl+-R2)(2fl^+21fli_R2-|-2-R^ 



Spherical scatterers 

Significant simplification is possible whenever A, B 
have spherical symmetry. First, the T-matrices are diag- 
onal in angular momentum basis and so may be expressed 
as 

2iLU . 



271 



2i5f 



where the normalization factor ^ follows from Eq (|47)) . 
A second consequence is that rotation around a (which 
from now on we take as coinciding with z axis) is a sym- 
metry of the whole system. The determinant therefore 
factors as a product of terms corresponding to different 
values of the azimuthal number m. The energy turns into 
a sum of the corresponding terms 



E ■ 



E 



duj 
2^ 



The matrices {k[)?^}'^,^^^^ defined for each m <^ "L 

(actually icC"") = are infinite dimensional but 

may be approximated in numerical calculations by finite 
matrices corresponding to < some Iq. The operator 
K'^'^^ may be written explicitly as 



\ogAei{l- K^"'\iuj)). 



(rn) 



(m)^(A) (m) (S) 



t, = i(-l)^(e2'*^-l). 



g'Cl (*^) = (-1)"+' VW+Wh + T) (60) 

k h l\( h hi 
/ I m -m 



Note that both g'f^i^[iLo) and tjiiuj) are real. 



Dirichlet Spheres 

The simplest example for which the above may be ap- 
plied is the interaction of two hard (Dirichlet) spheres. 
The T-matrix elements are well known in this case, and 



are given by tj(i^) — (— 1)^^""^ 



jj("-R) 



>{c.R) 



which translates to 



tj{iuj) 



(61) 



(59) on the imaginary frequency line. Here R is the sphere's 

radius and /, , 1 is a modified Bessel function of the first 

J I 2 

kind. 

In the special case where the two spheres have equal 
radii Ri = R2 an extra simplification occurs. One can 
then write K — K^, Kw — gii'ti', which imply logdet(l — 
K) = logdet(l + A')-|-logdet(l-^). The numerical cal- 
culation of the two determinants det(l±iir) is then some- 
what easier then direct calculation of det(l — K). More- 
over comparison with section IIVI shows that the two de- 
terminants det(l± if) (actually with i^ij' = (— l)™(7n'^i') 
correspond to the Casimir interaction energies Ejj^n of 
a single hard sphere and a Dirichlet/Neumann mirror a 
distance a/2 away. The symmetric two hard sphere sys- 
tem then has the energy Es = Ed + En- (One may also 
understand this in terms of decomposition into even and 
odd modes.) 

We have done the calculation including partial waves 



of I < Iq for different values of and considered the Iq 
dependence of the results as a test for convergence. Most 
calculations included modes up to Iq = 10, but for small 
values of sphere separation a, we used larger Iq even up 
to ^0 = 72 for a/R = 2.1. Since we expected the error 
to behave roughly as Ec — E{Iq) ~ 0[e~'^^°) we tried to 
fit the results with this assumed asymptotics. The num- 
bers suggests that in both the sphere-sphere and in the 
sphere-plate cases we have c ~ 2 log(l -I- d/R) where d is 
the distance between the two objects (i.e. d ^ a — 2R for 
Es and d = {a - 2R)/2 for Ed, Em). Table (p)) below 
shows the value of the constant c for identical spheres as 
a function of their distance as well as the value of at 
which the error dropped to within 1% of the exact re- 
sult. It should be remarked that the estimate for c is a 
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bit crude since our numerics is consistent with c being 
a slowly growing function of Iq (which might be due to 
sub- leading asymptotics). By matching our results with 
the assumed asymptotics, one can obtain a corrected es- 
timate for Ec- Comparison of this estimate with results 
obtained by increasing Iq gave good agreement. 



a/R 


c 


L{1%) 


2 





oo 


2.1 


0.18 


31 


2.2 


0.34 


16 


2.35 


0.57 


9 


2.5 


0.78 


7 


2.75 


1.06 


5 


3 


1.33 


3-4 


3.5 


1.78 


2-3 


4 


2.14 


2 


5 


2.75 


1 


7 


3.44 


1 



(62) 



The following table ([63]) and Fig[3] show the re- 
sults for the Casimir energy itself (measured in 
units of ^). Ed denotes the energy of Dirichlet- 
mirror+(Dirichlet)sphere system, En denotes the en- 
ergy of Neumann- mirror-|~(Dirichlet)sphere system, Es 
denotes the energy of the symmetric two hard sphere 
configuration (having Es — Eo + Em)- The result for 
Eb are in perfect agreement with a similar calculation 
done in Ref. [H. 



a/R 


Ed 


En 


Es 


2.1 


-8.75 


7.66 


-1.0939 


2.2 


-2.2129 


1.9382 


-0.27477 


2.35 


-0.739 


0.6488 


-0.0902822 


2.5 


-0.3688 


0.3245 


-0.0443005 


2.75 


-0.1679 


0.1483 


-0.0195891 


3 


-0.09703 


0.08613 


-0.0108937 


3.5 


-0.044981 


0.0403034 


-0.00467768 


4 


-0.0261973 


0.0236767 


-0.00252067 


5 


-0.0123048 


0.0112853 


-0.00101948 


7 


-0.00477708 


0.00447243 


-0.000304649 


10 


-0.00199796 


0.00190445 


-0.0000935083 


13 


-0.0011022 


0.00106165 


-0.0000405423 


16 


-0.000700129 


0.000678957 


-0.000021172 



(63) 



It may be remarked that Ed , E^ correspond to sphere- 
mirror distance, which is half the sphere-sphere distance 
in the corresponding calculation of Es- This fact is re- 
sponsible among other things to slower convergence in 
the calculation of Ed^En and hence to a smaller num- 
ber of calculated significant digits compared with Es- 




0.2 0.4 0.6 0.8 



FIG. 3: The calculated Casimir energy of: (a) Two Dirich- 
let spheres of radius R at distance a between their centers. 
(b,c) A Dirichlet sphere of radius R whose center is at a 
distance a/2 from a Dirichlet/Neumann mirror. The graphs 
show E/Eo as a function of {a/R — where Eo is the 

large distance asymptotic expression of it. Specifically: (a) 

En = --t-t-Kt^, (b,c) EF'^ = =F o , ^^T,vi ■ At short dis- 

4 

tances E/Eq approach the PFA prediction (a)^ ~ 0.27, 

(b)f^ ~ 0-54, (c)^ ~ 0.47. The black curve shows the 
calculated exact result for — > oo. We extrapolated it to 
a — 2R and o = cx3 using the known asymptotics. The col- 
ored graphs show the result of including partial waves oil <lo 
where: lo = Q (red), lo = 1 (sky blue), Iq — 2 (green), lo = 10 
(blue). 



We would like to mention two points regarding the 
actual implementation of the numerical calculation. (Our 
earlier numerical attempts failed because we were not 
fully aware of these points.) 

The expressions of giii[iuj),ti{iuj) may attain at small 
w's very large/small values respectively in such a way 
that only their product remain finite. At large w's 
similar phenomena occur with ti{iuj) large and gii>(iuj) 
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small. Thus to avoid computer overflow it is much 
better to "renormalize" these two quantities redefining 
gw = zizi'9w,ii = ti/z^ with zi - (i?a;)'+i/2e«". 

A second important point is that one should make sure 
that the computer program doing the calculation does 
not use the expansion of Ii^i/2{x) in terms of elemen- 
tary functions. In MATHEMATICA (which we used) 
this expansion is an automatic default whenever the in- 
dex of the Bessel function is half integer. However this 
expansion is known to be numerically unstable (except 
for very small I) and using it would lead to errors. 

The general formula works well for i?i 7^ i?2- For ex- 
ample taking i?i — Rq,R2 = 2Rq and measuring E in 
units of we found 



Interaction energy between a sphere of 
radius Rq and a sphere of radius 2Rq 



a/Ro 


E 


3.1 


-1.4554 


3.2 


-0.367535 


3.3 


-0.164591 


3.4 


-0.0931057 


3.5 


-0.0598295 


3.67 


-0.0334525 


3.83 


-0.021821 


4 


-0.01501511 


5 


-0.00362365 


6 


-0.00151965948 


8 


-0.00047970126 


10 


-0.00021536976316 


14 


-0.0000696380241 


18 


-0.00003103693506 


22 


-0.0000164921322 



(64) 



given by the following table ([65]) and Fig. [4l 



e 


771 

h 


64 


-0.003092 


100 


-0.003927 


900 


-0.00829 


lO'^ 


-0.0084835 


iU 


-U.Uiio4 


10^ 


-0.01364 


10^ 


-0.01447 


10^ 


-0.01483 


108 


-0.01495 


00 


-0.01501511 



(65) 




FIG. 4: The calculated Casimir energy of two (scalar) Dielec- 
tric spheres of radii R2 = 2Ri at centers distance a = 4Ri 
depicted as a function of their dielectric constant ei — £2- 
The energy E was normalized by the Dirichlet spheres result 
Eo so that at e — > 00 we obtain E/Eq = 1. 

The calculation may easily be repeated for any given 

IX. ELECTROMAGNETIC FIELD 

To extend the ideas of the previous section from the 
scalar to the EM case one needs to present the EM prop- 
agator in a form analogous to eqs. l(54|) - ((56)) . The re- 
quired representation of the EM propagator (derived in 
appendix]^ is 



^0 = \ {jma)B)Cj,na;j'ni'a'{{j'm'a' 



(66) 



Dielectric Spheres 



where a, a' can take the two values 0, 1 corresponding to 
the TE (magnetic multipole) or TM (electric multipole) 
modes respectively. The C coefficients are given by 



The formula also works well for finite dielectric con- 
stant. For example the numerical results for i?i = 
Ro', R2 — 2_Ro;a = 4i?o as a function of ei = £2 are 



C 



j7na;j'm' a' 



(67) 
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where may be defined in terms of vectorial spherical 
harmonics as 



(0) 



(1) - . J+1 



(68) 
(69) 



These functions satisfy iYj-^ —fx Yj^ , iY^^ = fx Y-^^ . 

After Wick rotating, we obtain Cjma;j'm'a'{i^) — 
p -j+a-a ^gj^^.j,^,^, where the coefficients 



In the particular case of a perfectly conducting sphere 
of radius R one has: 



'-TM 



(icj) 



X — LOR 



(74) 
(75) 



9jma;j'm' a' (i^) — \J u)a C^^) 

X E.'™- ^.'+1 (^«)^/"™" (a) / dm^n,' iYl,t ■ yS') 

are real. 

The integrals / dr2Y,3„3 (f/^"^* • Y^^J^) appearing in 
(|67l70p can be expressed explicitly in terms of 3j-symbols 
as follows. For a = a' it is given by 



h ( jl + 1) +32 (J2 + 1 ) - j3 (33 + 1 ) 

2Vii(3i+i)i2(i2+i) 
X i ] ( 

I M 



(2ji+l)(2j2 + l)(2j3 + l) 
47r 



(71) 



31 J2 J3 

-mi 7712 rris 



(which vanishes unless ji +j2 + js = 0mod2.) For a ^ a' 
the integral is nonzero only provided ji + j2 + j3 = lmod2 
in which case it is given by 




( 1)'"' ^/ (23l+l)(2j2 + l)(233 + l) ^J2) 
jl h j3 \ 

-TOi 1712 rns J 

Ea ([7T|) may be derived from Eq (|55)) by using the identity 
+ l)y,m ~ LYjrn (where L is the angular momen- 
tum operator) and integration by parts. The relation ()72p 
was found with the help of Eq. (18) in Ref. (35l |. 

Spherical scatterers 

Assuming spherically symmetric scatterers, one may 
define phase shifts S^rp^(uj), S^pj{u!) (by parity these two 
channels do not mix) . Similarly to the scalar case we use 
the notation t^o. = i (-l)J+«(e2'4°'(^-) - 1). 

Choosing the z-axis along d the operator K = TGTG 
splits to independent blocks corresponding to the 

values of the azimuthal number m. In a given block the 
g-matrix elements become 



(m) 
9ja;j'a' 



(73) 



j'm 



The matrix K'^"^\iijj) is then written explicitly as: 



K 



jcx.\j a. V / joi Jja;j a j a. ^ j ol \j a. 



(A)hyi) AB) (m) 



Using this we numerically calculated the electromagnetic 
Casimir energy for a pair of conducting spheres at dis- 
tance o between their centers. As in the scalar case writ- 
ing K = and considering det(l ± K) separately al- 
lowed us to also find the interaction energies Ee,Em of 
a sphere near a conducting/infinitely permeable mirror 
placed a/2 from its center. The two spheres energy is 
then the sum Eg — Ee + Em- Most of the calculations 
where done by including modes having j < 10, however 
for the shortest distances a = 2.35, 2.2, 2.1 were conver- 
gence is slower, we extended the retained modes up to 
j = 20,40,60 respectively. The results are shown in the 
following table (I76p (written in units where R = I) and 
in Fig. [5] 



a 




E,n 


Es 


2.1 


-16.15 




14.5 




-1.662 






2.2 


-3.82 




3.48 




-0.337635 




2.35 


-1.157 




1.073 




-8.356 


10 


-2 


2.5 


-0.53 




0.50 




-3.18 - 


10- 


-2 


2.75 


-0.211 




0.201 




-9.595 


10 


-3 


3 


-0.1074 




0.1036 




-3.787 


10 


-3 


3.5 


-3.97-10- 


2 


3.88-10- 


2 


-8.917 


10 


-4 


4 


-1.89-10- 


2 


1.86-10- 


2 


-2.864 


10 


-4 


5 


-6.24- 10- 


3 


6.19-10- 


3 


-4.887 


10 


-5 


7 


-1.38-10- 


3 


1.37-10- 


3 


-3.965 


10 


-6 


10 


-3.06 - 10- 


4 


3.06-10- 


4 


-3.032 


10 


-7 


13 


-1.04- 10- 


4 


1.04-10- 


4 


-4.703 


10 


-8 


16 


-4.47- 10- 


5 


4.47-10- 


5 


-1.085 


10 


-8 



(76) 



The numerical results seem to converge as jo ^ oo at 
roughly an exponential rate. The graph in Fig. [6] shows 
how the speed of convergence depends on the distance 
between the bodies. It is interesting to note that the 
results obtained in section IVIIII for the scalar case give 
almost the same graph. Also, one can easily check that 
the results for Eg are basically the same as the ones ob- 
tained in [2lj taking into account we chose to normalize 
the energy in comparison with the large distance asymp- 
totic expression for the energy. 
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0.2 0.4 0.6 0.8 



FIG. 5: EM Casimir energy of two conducting spheres of 
radius R at distance a between their centers. The graphs 
show E/Eq as a function of (a/i?— 1)^^ where Eo is the large 



143ii'^ 
' 16ir(a-2i{)^a 



distance asymptotic expression of it. Eg 
The black curve shows the calculated exact result for jo ^ oo. 
We extrapolated it to a = 2R and a = oo using the known 
asymptotics. The colored graphs show the result of including 
partial waves of j < jo where: jo = 1 (red), jo = 2 (sky blue), 
jo = 4 (green), jo = 10 (blue). 
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3 

2 
1 

d/R 



10 



12 



14 



FIG. 6: Including partial waves of j < jo results in error 
behaving roughly as e"^-'". The graph shows the constant 
c as a function of the separation distance d. The blue dots 
correspond to our results for two conducting spheres (where 
d = a — 2R) and the red dots to conducting-sphere + con- 
ducting plate (where d = a/2 — R). At small distances both 
cases give c ~ 1.1 d / R. 



APPENDIX A: PROOF OF THE G REEN' S 
FUNCTION EXPANSIONS EQNS (|54|66|) 



Scalar case 



Suppose R = a -\- •f then obviously e 
Inserting the well known expansion 

e*-^ = 47r^*'y,:„(fc)11,„(f)j;(fcr 

We get 



ik-R ^ik-a ^ik-r 



e e 



Multiplying both sides by Yira{k) and integrating J dftk 
we find 



ji{kR)Yi^{R) = 
X ika)Yl„,r^" {a)jv ikr')Yl,,r^' (r') 



(A.l) 



Concentrating on the case R,a > r' , it makes sense to 
separate the ingoing and outgoing parts in the last equa- 
tion. This amounts to replacing the bessel functions 
ji{kR), ji{ka) by hankel functions hi{kR), hi{ka) of the 
first or second type corresponding to outgoing or ingoing 
waves. Since this argument may seem as hand-waving, 
we will return and elaborate on it more at the end of the 
proof. Equating the outgoing parts we have: 



h[^\kR)Yi^{R) = 



(A.2) 



h\:>{ka)Yi.,.^,.{a)ji.{kr')Yi,^.{f') 



It is well known that for R > r the free propagator 
may be expanded as 



1 



1 



^ik\R-r\ 



Substituting here Eq. (|A.2p we finally get 



kJ23iikr)h['\kR)Y;^{f)YUR) 



G(f,a 



(A.3) 



,(1) 



U dnY,„,Y*^.Y*,^„) 



xji iu;r)ji, {u;r')h\t> (wa)^',™ ir)Yi'm' (r')l;"m" (a) 



Which is exactly Eq. (f54|) . 

Let us now return to the derivation of Eq. (|A.2|) 
from Eq. (|A.1|1 . We first note that the function 



h[^^ {kR)Yigmo (R) with R — a + f being a solution of the 
free wave equation may be expanded around r = in the 
form 



(1) ~(i) 



To be more 



for some (a dependent) constants c| 
precise {kR)Yi„ma {R) is a solution only for f ^ —a 
(i.e. R 0) therefore one has two separate expansions: 
one for r < a and another for r > a. We concentrate on 
the latter. 

Since h'l^\kR)Yigmg{R) is a purely outgoing wave it 
is clear that the expansion in terms of r must also con- 
tain only outgoing waves i.e. = 0. This claim is 
based on "physical intuition" . A more rigorous mathe- 
matical argument may be constructed by considering first 
pure imaginary k — iq with q > 0. One then note that 
h^j^ [iqR) is exponentially decreasing as i? — > oo which 
imply that the same must hold for the r.h.s. Since the 
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Yim 's are linearly independent this require all the 's 
to vanish. 

A similar expansion obviously exists also for /i'^' : 
h^JikR)Yi„„,„{R) ^Y.^?r!.hf\kr)YUf) 
Summing the two expansions we have 

ji,{kR)Yi,^,(R) ^ \ {h\l\kR) + hl^\kR)) = 

However such an expansion is clearly unique. Therefore 
it must be the same as the expansion in Eq. (|A.1[) . Com- 
paring the two (and using ji = \ {h^^'> + /i*-^-*)) we deduce 



Irn Irn 



47r (j dmoraoYLYv'r 



which proves Eq. (|A.2[) . 



xji„{ka)Yi„„i„{a) 



Here Yjim are vectorial spherical harmonics, Yim are the 
usual scalar spherical harmonics, and j/, hi are spherical 
Bessel and Hankel functions. The coefficients 




are found to be expressed as the following integral of 
spherical functions: 



I I' I" \ _ 
m m! m" j 

^7rJdn{Y,i„,-Yp,„,)Y*^„ 



(A., 



The radiation gauge propagator Vq is given by the 
transverse part of €So. In Eq (jA.6[) each j,m corre- 
spond to three different spherical function \jlm) (having 
I ~ j ~ 1, j, J + 1). These may be decomposed in terms of 
the TE and TM modes and a nonphysical longitudinal 
mode. 

\TE) = \jjm) 



\TM) 



2j + l 



2j + l 



+ l,m) 



The electromagnetic case 



To derive the EM expansion ([66)1 we similarly start by 
using the identity 



^^Y.^'3i{kr)Y;i„,{k)®Y,i„,{f) 



(A.4) 



Repeating the same steps as for the scalar we then find 
that 

Go = X 1 

Airr 

may be expanded as 

= \{]lm)B)C.ji^;ri'^'{{j'l'm')A\ (A.5) 

where 



2uj 



\{i^m)AM) = \j —il{^rA,B)Yjlm{fA,B) (A. 6) 

are the free vectorial spherical wave functions centered 
at Pa , Pb ■ The C coefficients may be written as 



C 



jlm:j' I' m' 



ITT 

2uj Z-//",m" 



(A.7) 



m m m 



h\,/{uja)Yi»rn"{a) 



2j + l 



^• + 1 

2j + 1 



+ l,m) 



To obtain the required expansion of the radiation 
gauge propagator T>q we need to rewrite Ea (|A.5p in terms 
of these three modes and drop the parts containing the 
longitudinal mode. This can be done quite straightfor- 
wardly leading to the results (|66ll69p . 



APPENDIX B: ANALYTICAL PROPERTIES OF 

THE TaGoTsGo-OPERATOR: 

Having established the form ([TT|) for the energy, we 
turn here to discuss the properties of this expression. 
The main aim of this appendix is to rigorously show that 
the object logdet(l — TaGoabTbGoba) is well defined 
and finite. The main mathematical notions and theorems 
which we use here, are briefly reviewed in appendix ICl 

As already remarked in the introduction it is well 
known that det(l — M) is well defined whenever M is 
a trace class(t.c.) operator (definition [C3}. We would 
like to show that for a large class of situations (including 
a pair of disjoint finite bodies A, _B, separated by a finite 
distance) the operator TaGqabTbGqba ■ Ha Ha is 
trace class in the continuum limit , and so prove that 
indeed the expression (|lip is finite and well defined. 
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Indeed, by theorem IC. 51 the mere fact that Go{x,y) is 
a smooth function for x y is sufficient to guarantee 
that for any pair of compact volumes A, B gM.^ at finite 
mutual distance the operator Gqab is trace class. To de- 
duce that TaCqj^bTbGoba is trace class (and by similar 
argument also 1 — GqJTa appearing in Eq. (|^ ) it is 
then enough (proposition lC.6]) to make sure TA,B{i^) are 
bounded (definition IC.2p . 

In the context of dielectric interaction, it is particularly 
easy to show that T{iuj) is bounded. In physical systems 
at equilibrium, it follows from causality properties of the 
dielectric function [3l[, that x(ia;, x) > 0. We then have 
the following 

Lemma B.l. For x(iLU,x) > , the T operators are 
positive and hounded. 

Proof: Since GcX > (definition IC.3|) one may write 

2 

T — ^Jx ^_^j2.'^Q^^ \/x from which it is seen that T > 

and that in the operator norm ||r|| < cij^||x|| □. 

In fact, this holds also for nonlocal x as long as /(x) ^ 

x(*'^j 2:, x')/(a:;')da:;' is a bounded positive operator 
Ha — > Ha- In the context of more general type of inter- 
actions which may not be positive, one needs to use some 
assumption on the stability of the system to guarantee 
that TiiiS) is bounded. We do not elaborate on this here. 

An alternative approach to proving the trace class 
property of TaGqabTbGqba is based on the notion of a 
Hilbert-Schmidt operator (definition IC.7[ also denoted 
H.S.). Here the frequently used strategy in operator 
analysis is to use the following fact: HUG H.S. and 
V £ H.S., then UV G t.c. The advantage of this ap- 
proach is that it is very easy to check if an operator 
is Hilbert Schmidt. Since the Hilbert Schmidt norm is 
ll^^lllfs ~ Tt{A^A), one may evaluate it directly, (e.g. 
by computing J \A{x,x')\'^). 

Theorem B.2. For any two bodies A, B such that 
X4xB^^^y|'^o(a;,y)P < oo, TaGoabTbGqba is trace 
class. 

Proof. First we show that TaGqab and TbGqba are 
Hilbert Schmidt operators. This can be verified in the fol- 
lowing way. We have just seen that Ta,Tb are bounded 
operators. Now note that Gqab is Hilbert-Schmidt, 
since, 

WGqabWh.s ^ dxdy\GoABix,y)\'^, (B.l) 

Jaxb 

which is finite under the condition above. Now the 
inequality \\TaGoab\\h.s < \\Ta\\\\Goab\\h.s implies 
that TaGqab is Hilbert Schmidt. Finally using U,V £ 
H.S. ^UV £ t.c. we see that TaGqabTbGqba e t.c. □ 

Corollary B.3. For any finite bodies A, B, such that 
distance{A, B) > 0, and any Green's function which is 
finite away from the diagonal, TaGqTbGq £ t.c. 




FIG. 7: The support of the function /a 

Example B.4. For the scalar field discussed above, 
GQ(x,y) — 4jr|l_J| 7 condition is satisfied. In the 
same way it is satisfied for the electromagnetic field (one 
has to take into account also matrix indices but these dis- 
crete indices do not change finiteness of the integrals ) 

Remark B.5. The uj integration in Eq. (jll[) . is con- 
vergent. To see this note that Go decays exponentially 
with oj therefore, \\Gq\\h.s. decays exponentially, also the 
\\T\ \ 's do not grow more then quadratically in uj. 

In the EM case one may also worry due to the factor 
^ appearing inVoijix,y) = {6,^^ - ^V'f^V''f^)Go{x,y), 
about convergence for to ^ 0. This factor however gets 
cancelled since \\T\ \ < a;^||x|| as shown in lemma VB.U 

One may also show that Gqab are t.c. themselves by 
using H.S. properties. The bodies are assumed not to 
touch, thus we can choose a (compactly supported 
and infinitely smooth) function such that PAfA — 
Pa, and PBfA = where Pa, Pb are the projections on 
L'^{A),L'^{B) (i.e. /^(x) = 1 for x £ A, and it then 
smoothly goes to 0, before reaching body B see Fig [7]) . 

Writing: 

GoAB = L1L2 (B.2) 

il = Paj^^t^ ; L2 = (p2 + oJ^YfAG^PB, 

we see that if 4a > rf, 

II^iIIh.s. = '^^{Paj^^^){Paj^^^V = (B.3) 
Vol{A)jd<'p\j^^,:^\<^ 

and so Li is Hilbert Schmidt. Next, we check that L2 £ 
H.S.. To see this last point, note that 

<x\L2\x' >=< x\{p' + uj'')'^ fAGoPBW >= (B.4) 
(-A, -I- u;^rfA{x)Go{x - x')Pb{x') 

Since Go (a; — x') is smooth away from x — x' , where 
the expression is anyway zero because fAPB = 0, and 
since {x\L2\x') has compact support (for integer a) we 
see that 11^2111/ s = / da;da;'|L2p < 00. Thus, Gqab can 
be written as a product of two H.S. operators, and as 
such is trace class. 
Finally, we have that 
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Theorem B.6. ^Eigenvalues of TGTGj For x > 

0, all eigenvalues A of the (compact) operator 
TaGqabTbGqba appearing in pT|) satisfy 1 > A > 0. 

Proof: Wc will use repeatedly that for bounded opera- 
tors X, Y the nonzero eigenvalues of XY and YX are the 
same. Note first that GojX > (as operators) implies 

spec(xGo) \ {0} = spec{y/G~oxVG^>) \ {0} C [0, oo). 

Writing TaGo = 1 — as an operator on L^(R'^) 

it is then clear that its spectrum lies in [0,1). The same 
conclusion then applies to the operator ^/GoTa^/Go but 
since it is hermitian one concludes also 1 1 VGoTli VGo| | < 
1 from which it follows WVG^TaGqTbVG^W < 1 and 
hence A < 1. Similarly ^/U^T^V^ > imply A > □ 



Theorem C.5. Consider an operator A : L'^{Di) 
L^{D2) where Di,D2 are some domains in R" 
which is given explicitly as an integral A'ip{x) = 
K{x, y)ip{y)dy. A sufficient condition for A to be 
trace class is that Di,D2 are compact and K{x,y) is 
smooth in a neighborhood of Di x D2 ■ 

Proposition C.6. // A is trace class and B 
bounded then AB and BA are also trace class and 
Tr{\AB\),Tr{\BA\)<\\B\\Tr{\A\). 

Definition C.7. M is a Hilbert Schmidt operator (de- 
noted M e H.S. or M e J2) if I \M\ ^ ee TrMU/ < 00 

In particular we mention that the product of two 
Hilbert Schmidt operators always give a trace class op- 
erator. 



APPENDIX C: SOME PROPERTIES OF 
(INFINITE DIMENSIONAL) OPERATORS 

Here we recall some mathematical notions that we have 
used in describing the trace class properties of Eq. pT|) . 

Definition C.l. For an operator B : H H , 

the operator norm of \\B\\ is defined as \\B\\ = 

\<^\B\H,>\ 
sup^g/i-,^,^0 <V)|v.> 

Definition C.2. An operator B is bounded if \ \B\\ < 00 

Definition C.3. An operator A : H ^ H is called a 
positive operator (denoted A > 0) iff {ip\A\iJ;) > for 
every ijj H . 

This implies that A is hermitian and its spectrum non- 
negative, li A : H ^ H is a. positive operator then there 
exist a unique positive operator B : H ^ H satisfying 
A = B^. B is called the square root of A and denoted 

Va. 

Definition C.4. An operator A : Hi ^ H2 is called 
trace class (and denoted A E t.c. or A £ Ji) iff 

^ ll^V'nll < 00 where {'ipn}'^=i is some orthonormal ba- 
sis of Hi . It can be shown that this condition does not 
depend on the choice of the orthonormal basis. (Note that 
the definition makes sense even when Hi ^ H2-) 

li A : H ^ H is trace class then for any orthonormal 
basis {'0n}^i of H the sum ^(^/;„|^|i/'n) converges to 
the same (finite) value which is denoted tr{A) and called 
the trace of A. One then also have tr{A) = where 
{A„} are the eigenvalues of A (Lidski's theorem) 

If A : H ^ H is trace class then the determinant 
det(l + A) may also be rigorously defined and one has 
det(l + A) = n(l + A„). 

The following theorem may be proved using the well 
known fact that the Fourier coefficients of a smooth 
K{x,y) decay faster then any power. (Note that these 
coefficient also serve as the matrix elements with respect 
to Fourier basis of the operator defined hy K.) 
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